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Abstract. The mass transport in soft-sphere mixtures of small and big particles 
as well as in the disordered Lorentz gas (LG) model is studied using molecular 
dynamics (MD) computer simulations. The soft-sphere mixture shows anoma- 
lous small-particle diffusion signifying a localization transition separate from the 
big-particle glass transition. Switching off small-particle excluded volume con- 
straints slows down the small-particle dynamics, as indicated by incoherent in- 
termediate scattering functions. A comparison of logarithmic time derivatives of 
the mean-squared displacements reveals qualitative similarities between the lo- 
calization transition in the soft-sphere mixture and its counterpart in the LG. 
Nevertheless, qualitative differences emphasize the need for further research elu- 
cidating the connection between both models. 

1 Introduction 

Transport in heterogeneous disordered materials and in porous media arises in many diverse 
subjects of science and engineering such as water resource management, oil recovery, the physics 
of crowded biological systems, glass technology or even the geophysical understanding of the 
eruption of volcanos Examples are fast ion transport in alkali-silicate melts and other 
industrially relevant glass formers , protein motion on cell membranes [3, |3| and within the 
living cell among others. Common to all these materials is that they consist of at least 
two components, one of which (the "fast component") is responsible for the observed transport 
phenomena, while the other forms the heterogeneous matrix and is either completely frozen as, 
e.g., in porous media, or relaxes orders of magnitude slower, as in glass- forming ion conductors. 

A well-known reference point for diffusion in random environments is the Lorentz gas 
(LG) [6, 7]: considering a single particle diffusing through a set of fixed obstacles, one essentially 
arrives at an off-lattice model for the transport through porous media or ion-conducting sys- 
tems. The LG exhibits a localization transition where diffusion of the tracer ceases altogether if 
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Fig. 1. Incoherent intermediate 
scattering function Fs{q,t) for the 
small/tracer particles. Solid and 
dashed lines correspond respectively 
to the BSSM with and without 
interactions among the small par- 
ticles, while dotted lines represent 
the LG model. For the BSSM, 
the wavenumbers are in the range 
0.74 < qcT < 17.5, as indicated (with 
a the diameter of the big particles). 
The wavenumbers in brackets corre- 
spond to the LG. The number density 
of the BSSM is p = 3.26. For the LG, 
the packing fraction of the obstacles 
is n* = 0.75. 

the density of scatterers exceeds a certain value; this transition is understood as a dynamic crit- 
ical phenomenon connected with the percolation of void space in the heterogeneous background 
medium 

Obviously, several aspects of ion-conducting systems are described in an oversimplified man- 
ner by the LG model. In a real ion conductor, the frozen environment is not formed by fixed 
obstacles that are randomly distributed. Instead, the "obstacle particles" exhibit structural 
correlations (at least on short length scales) and they are subject to thermal motion. Moreover, 
in most ion-conducting systems the mobile ions cannot be considered as non-interacting tracer 
particles and the correlations between ions have to be taken into account. Realistic modeling for 
sodium silica melts [Tol - [T2l | reveals that the silica network can constitute such a quasi-arrested 
array of correlated obstacles through which sodium ions meander on preferential diffusion path- 
ways. The chemical properties of the silica melt to form a tetrahedral network appears not to 
be essential: a similar decoupling of diffusive transport is observed in dense size-asymmetric 
Yukawa melts [l^ and binary soft-sphere mixtures fl3 . Mode-coupling theory has been used 
to predict it for size-disparate hard-sphere systems |l5l . IT6| . As shown recently by molecular 
dynamics (MD) computer simulations [l7l |. such systems may exhibit a glass transition where 
the fast component remains mobile even for high-density states at which the matrix does not 
show relaxation over the entire simulation time window. At these high densities, evidence was 
found that the small particles show anomalous diffusion and approach a localization transi- 
tion that bears resemblance to the LG transition. Similar findings hold for an entirely frozen 
matrix, where non-trivial predictions of a mode-coupling theory developed by Krakoviack [l^ 
have been tested recently [l^, However, it is an open question to what extent the small- 
particle localization dynamics in binary soft-sphere mixtures can be understood as a dynamic 
critical phenomenon that falls into the same universality class as the localization transition of 
the LG [S, i9|] and other continuum percolation models ^21i] . In particular, the question arises 
how the localization of the small particles is affected by interactions between the mobile tracers. 

In this contribution, we present first steps to address the latter issues. We compare on a 
qualitative level the anomalous diffusion of a binary soft-sphere mixture with that of a LG. The 
effect of collective interactions among the small particles is investigated by switching off their 
interactions. 



2 Results 

We performed MD simulations of an equimolar binary mixture of purely repulsive soft spheres 
(BSSM) with a size ratio of 0.35. Diameters are chosen additively, and nonadditive energetic 
interactions further decouple the species. Temperature is unity, and all masses are equal. For 
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Fig. 2. Effective exponent, as obtained from the 
logarithmic derivative 7(t) of the mobile parti- 
cles for the binary mixture with and without in- 
teractions between the small particles as well as 
for the LG. Results are shown for different den- 
sities: the three sets marked with number 1 cor- 
respond to p = 2.296 and n* — 0.40, those with 
number 2 to p = 3.257 and n* = 0.75 and those 
with number 3 to p = 4.215 and n* — 0.82 (here, 
p and n* are the densities of the binary mixtures 
and the LG, respectively). 



details, see Ref. [T3- MD simulations of the LG were carried out using an event-driven algo- 
rithm and randomly distributed, possibly overlapping hard sphere obstacles. The dimensionless 
obstacle density n* is the only control parameter of the model; details can be found in Ref. |^. 

Figure [T] shows the time dependence of the incoherent intermediate scattering function, 
Fs{q,t), of the small particles for different values of the wavenumber q, as indicated. For the 
BSSM, at a number density p = 3.26, the big particles form a glass and the transport of the small 
particles is characterized by anomalous subdiffusive motion crossing over to ordinary diffusion 
at long times. Two kinds of BSSM systems are considered in Fig. [1] the solid lines correspond 
to the system with interacting small particles, while the dashed lines correspond to a system 
where the small-small interaction potential is set to zero, keeping all other parameters fixed. 
Structure and dynamics of the big particles as measured through the partial static structure 
factor are unaffected by this jlTj . Surprisingly, at fixed q the curves for the fully interacting 
case decay more rapidly than those for the non-interacting case. In Ref. 17, we have proposed 
that this difference arises because small-particle interactions result in a more directional and 
thus more effective diffusive motion. 

Also shown in Fig.[l]are results for the LG at a packing fraction of the obstacles n* = 0.75. 
At this value of n*, the decay of Fs{q,t) takes place on a similar time scale provided that 
one reduces the wavenumbers by about a factor of 4 (as done in Fig. [T|) to match the length 
scales of the LG and the BSSM systems. The curves for the LG demonstrate the absence of 
a "glassy" plateau at intermediate times. At long times, they decay to a plateau of different 
origin: the presence of closed, finite pockets in the obstacle matrix leads to a localization of 
tracer particles. However, apart from that the decay of the LG correlators is similar to that for 
the BSSM systems. So the question arises whether the transport as seen for the non-interacting 
particles in the binary mixture is similar to that seen in the LG. 

To address this question, we compare the logarithmic derivative j{t) — d[log(Sr'^(i)]/d(logt) 
for the two binary mixtures and for the LG in Fig. [2] j(t) is calculated from the mean-squared 
displacements Sr'^{t) of the mobile particles in both systems. The quantity j(t) has the meaning 
of an effective exponent that crosses over from 7(0) = 2 (ballistic short-time motion) to 7(00) = 
1 for diffusive particles in the liquid. For localized particles, 7(00) = is expected. 

The selected obstacle densities for the LG in Fig. [51 n* = 0.40, 0.75 and 0.82, are below the 
critical density of the localization transition, n* — 0.837 Q. At criticality n* — n*, the exponent 
7(i) « 0.32 is approached at long times @, thus confirming the universal subdiffusive behavior 
expected from renormalization group arguments (refer to Ref. for a detailed discussion). 
Below n* , such anomalous diffusion is seen in a finite, intermediate time window and universal 
corrections to scaling lead to the observation of an effective exponent that is larger than the 
universal critical exponent. At n* = 0.82 for example, an effective exponent 7 w 0.4 is observed 
over about four orders of magnitude in time, see Fig. [21 
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At low densities (set of curves "1" in Fig. [2]), the binary mixtures with and without interac- 
tions between the small particles show almost identical behavior as the LG with respect to "f{t). 
Nevertheless, there are quantitative differences: At higher densities the effective exponent of 
the BSSM without interaction becomes rather small, the corresponding mean-square displace- 
ment remains of the order of the interparticle distance, and the cage effect interferes strongly 
with the regime of anomalous diffusion. Turning on the interaction, a window of subdifFusion 
emerges, yet the effective exponent drifts gradually from 0.5 to 0.6. In contrast, the LG displays 
clear subdiffusive behavior over several orders of magnitude in time as manifested in an almost 
constant 7(i), as mentioned above. 

Clearly, a simplistic mapping between the BSSM and the overlapping LG is not obvious 
and may not exist due to the highly simplified character of the LG. Thus, the investigation 
of LG variants with increasing complexity is desirable, which close the gap to the BSSM. For 
example, a better approximation of the matrix structure would be obtained by introducing 
correlated obstacles in the LG ("non-overlapping LG"). Another issue is raised by the hard 
interaction potentials in the LG. While a tracer with constant velocity between soft obstacles is 
easily mapped to the hard sphere model, the percolation transition is already smeared out for 
thermal tracers (or tracers in contact with a heat bath) due to the average over the Maxwell 
distribution. All these issues are the subject of forthcoming studies. 
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